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^ : ABSTRACT 

5h 

Thermal-gravitational instability (TGI) is studied in the protogalactic envi- 
^ ■ ronment. We extend our previous work, where we found that dense clumps first 

form out of hot background gas by thermal instability and later a small fraction 
>- ■ of them grow to virialized clouds of mass Mc > 6 x 10^ M© by gravitational infall 

and merging. But these clouds have large angular momentum, so they would 






•/^ ' be difficult, if not impossible, to further evolve into globular clusters. In this 

Q ■ paper, through three-dimensional hydrodynamic simulations in a uniformly ro- 

^ ■ fating frame, we explore if the Coriolis force due to rotation in protogalactic disk 

■^ ' regions can hinder binary merging and reduce angular momentum of the clouds 

Q^' formed. With rotation comparable to the Galactic rotation at the Solar circle, 

Q ■ the Coriolis force is smaller than the pressure force during the early thermal in- 

-<— > ' stability stage. So the properties of clumps formed by thermal instability are 

ci • not affected noticeably by rotation, except increased angular momentum. How- 

ever, during later stage the Coriolis force becomes dominant over the gravity, 
^ ■ and hence the further growth to gravitationally bound clouds by gravitational 

j^ ■ infall and merging is prohibited. Our results show that the Coriolis force effec- 

tively destroys the picture of cloud formation via TGI, rather than alleviate the 
problem of large angular momentum. 
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1. Introduction 

Thermal instability (TI) (Field 1965) is often invoked to explain a variety of physical 
phenomena in astrophysical environments: for instance, the multiple phases of interstellar 
gas {e.g., Field et al. 1969; McKee & Ostriker 1977), the formation of globular clusters 
{e.g., Fall & Rees 1985; Kang et al. 2000), cooling flows in clusters of galaxies {e.g., Nulsen 
1986), and the generation of turbulent flows in the interstellar medium {e.g., Vazquez- 
Semadeni et al. 2000; Kritsuk & Norman 2002). Although the basic concept of TI as a local 
instability is rather simple and robust, the realistic situation is often more complex, involving 
effects such as magnetic field, turbulence, gravity, and rotation in case of galactic disks. 

Recently the formation of structures via thermal-gravitational instability (TGI) in the 
protogalactic halo environment was studied using three-dimensional numerical simulations 
(Baek et al. 2005, hereafter Paper 1). The growth of density perturbations initially via TI 
and subsequently via gravitational infall and merging was followed up to 16 — 20 cooling times 
in a periodic cubic box with size L = 10 kpc. The simulations showed that clumps emerge 
first on scales smaller than the cooling length as a result of non-linear behavior of TI. Those 
clumps grow through compression by background pressure, as well as through gravitational 
infall. Later during the gravitational merging stage, some clumps become gravitationally 
bound, virialized clouds with mass M^. > 6 x IO^Mq and radius Re ~ 150 — 200pc. However, 
these massive clouds acquire angular momentum through tidal torque and merging and have 
a large spin parameter (A^) ~ 0.3. Hence removal of the angular momentum from theses 
clouds is critical, if they were to collapse further to form halo globular clusters such as in 
the model by Fall & Rees (1985). 

In this paper we study the effects of rotation in protogalactic disk regions on the for- 
mation of clouds via TGI. Uniform rotation was included in the same numerical simulations 
as in Paper 1 to model protogalactic disk environment. In the next section we describe 
our models and numerical method. The simulation results are presented in §111. Summary 
follows in 5IV. 



2. Numerical Simulations 

As in Paper 1, we consider a primordial gas of Th = 1.7 x lO^K, which corresponds to 
the canonical temperature of an isothermal sphere with circular velocity Vc = 220 km s~^, 
representing the hot phase of the gas in disk galaxies like the Milky Way. The fiducial value 
of the mean background density of hydrogen nuclei was chosen to be n/j = 0.1 cm~^. For 
the primordial gas with an assumed ratio of He/H number densities of 1/10, the gas mass 
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density is given by ph = (2.34 x 10"^'' g) Uh- With T/^ = 1.7 x lO^K and nu = 0.1 cm~^, 
the initial coohng time scale is tcooi = 2 x 10^ yrs. On the other hand, the free-fall time 
scale, or the gravitational time scale, is tgrav = 1.4 x 10^ yrs ~ 7 tcooi- Note that tcooi ex ^h^ 

— 1/2 

while tgrav ^ iT'h ■ ^'^ cooling, compared to gravitational processes, becomes relatively 
more important at higher densities and vice versa. The cooling length scale is defined as 
^cooi = Ch ■ tcooi = 4 kpc, where Ch = 198 km s~^ is the sound speed. 

To model the protogalactic disk environment, uniform rotation was incorporated by 
using a Cartesian coordinate system that rotates with angular velocity, Qz. We chose as 
a fiducial value for the angular speed ilo = 27 km s~^ kpc^^, which corresponds to the 
Galactic rotation at the Solar circle {Ro ~ 8.5 kpc) (Feast & Whitelock 1997). A case with 
lower the angular speed, fio/2, was also considered for comparison. In the rotating frame, 
an additional Coriolis force term, /,, = —2Qz x v, is added to the equation of motion. The 
simulation box was set to be cubic, with size L = 10 kpc = 2.5/cooi- The size was chosen 
to be large enough to accommodate a fair number of thermally unstable clouds of cooling 
length size, and so to obtain good statistics of cloud properties. Period boundary condition 
was imposed on the box, although it might not be the most natural choice for a disk-like 
geometry. We believe this particular choice of boundary condition would not affect the main 
conclusion regarding the role of rotation. Our periodic, cubic, simulation box represents a 
volume of disk region with significant rotation inside a protogalaxy. 

To mimic density perturbations existed on a wide range of length scales inside the 
protogalaxies, the initial density field was drawn from random Gaussian fiuctuations with 
predefined density power spectrum. The density power spectrum was assumed to be given by 
Pk oc k"' with n = (white noise). In Paper 1 we showed that the properties of clouds, once 
formed, do not depend on initial fluctuation spectrum, although their spatial distribution is 
sensitive to the spectrum. Without a priori knowledge on initial fluctuations, a spectrum 
of constant power over all scales was chosen. The amplitude of the density power spectrum 
was flxed by the condition 6rms = (^P^) / (p) = 0-2- The initial temperature was set to be 
uniform, and the initial velocity was set to be zero everywhere. 

The evolution of the gas from the initial perturbations was followed with 1) radiative 
cooling due to the primordial gas (Sutherland & Dopita 1993) down to T = 10"^ K, 2) 
background heating equal to the cooling of the initial, unperturbed background gas, and 
3) self-gravity. The hydrodynamics was solved using an Eulerian hydrodynamic code based 
on the total variation diminishing scheme (Ryu et al. 1993) on a grid with 1024^ cells (or 
512^ cells in the comparison run with Qo/2). Simulations started at t = and lasted up to 
tend = 16 — 20tcooi- Three simulations are presented in this paper, differing in angular speed. 
Model parameters are summarized in Table 1. 
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3. Results 

We begin our discussion by comparing the Coriolis force to pressure force and gravity 
for the flows associated with the formation of clouds. First, the ratio of the Coriolis force to 
pressure force can be estimated roughly as 

/e 2noV lOQoRc .,, 

rsj r\-/ I J^ J 

fp Ph/RcPh 3c/i 

where the typical flow speed is assumed to be similar to the sound speed, i.e., v ~ Ch, and 
the adiabatic index of gas is 7 = 5/3. Here Re is the radius of typical clouds. Conservatively 
with Re ~ 100 pc, the ratio is fc/fp ~ 0.05. So we expect that the effects of rotation are 
small, if not negligible, during the TI stage. The ratio of the Coriolis force to gravity can be 
written roughly as 

k ^ ^3^ (2) 

Again with v ~ Ch, and conservatively with Mc ~ 10® M0 and Re ~ 100 pc, the ratio 
becomes fc/fg ~ 25. Hence during the gravitational infall and merging stages, the Coriolis 
force is expected to play a dominant dynamical role in our simulations with rotation. 

The Coriolis force causes flows to be deflected at the right angle to the flow direction, 
hence it hinders spherical infall motions toward high density peaks, and generates circular 
motions in the plane perpendicular to the rotation axis. The Coriolis force also hinders 
the mergers of two clumps by drifting them to the opposite directions perpendicular to 
the line connecting two clouds, which leads to stretched worm-like structures. As a result, 
the formation of knot-like structures is suppressed and instead filamentary and sheet-like 
structures seem to appear. 

Figure 1 shows the density power spectrum at different times. In the figure the di- 
mensionless wavenumber is given as A; = L/\. The power spectrum is presented in a way 
that f Pkdk = (p^). During early stage t < 6tcooh the evolution of the density power spec- 
trum looks similar in two models S1024 and R1024. It is because with tcooi < ^grav, initially 
the power grows mostly due to TI, and the Coriolis force has negligible effects. During 
t ^ Stcooij gravity controls the growth and so the two models evolve differently. In model 
S1024 the power continues to grow due to gravitational infall and merging over all scales. In 
model R1024, on the other hand, the Coriolis force suppresses the growth, more evidently 
on smaller scales. The comparison of the power spectrum of models R1024 {Q = Qo) and 
RH512 {Q = O.5Q0) demonstrates the effects of different rotation speed. The figure shows 
that even with a half rotation speed the Coriolis force suppresses the growth of the power 
spectrum effectively. 
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To look at the shape and spatial distribution of formed structures, three-dimensional 
isodensity surfaces of a cubic region are shown in Figure 2. Initially clumps appear via TI 
in both models, as discussed in Paper 1 (see Figure 3 there). In model S1024, by t = Stcooi 
the clumps have developed into distinct clouds with roughly spherical shapes. By t = 12tcooi 
some clouds have grown to be massive enough to become gravitationally bound, and by 
t = 16tcooi they have grown more massive. Those gravitationally bound clouds have the 
central density higher than 1000 times the mean background density. In model R1024, on 
the other hand, the initial clumps do not grow to distinct clouds. Instead loosely connected 
filamentary and sheet-like structures appear with the central density lower than 100 times 
the mean background density. 

Although clouds are not distinctively defined in model R1024, we still identified clouds 
around high density peaks by the algorithm clumpf ind described in Paper 1, and calculated 
their properties. The first row of Figure 3 shows the number of identified clouds, Nc, as 
a function of their mass, Mc. The mass function of the two models is almost identical at 
t = 4tcooi- It is roughly Gaussian, since the initial density perturbations were drawn from a 
random Gaussian distribution. In model S1024, the mass function has evolved roughly into 
a log-normal distribution by t = Stcooh which is a signature of nonlinear structure formation. 
At later stage the mass function extends to higher mass with a power-law distribution, as 
more massive clouds develop through gravitational merging. In model R1024, however, the 
mass distribution remains roughly Gaussian at t = Stcooi, and later it develops into a form, 
which is not well defined. This is another indication that the development of nonlinear 
structures has been severely altered by the Coriolis force. 

The second row of Figure 3 shows the energy ratio of identified clouds, f3 = 2{Ek + Et)/\Eg\ 
Here E^ is the kinetic energy defined in the center of mass of a given cloud, E^ is the ther- 
mal energy, and Eq is the gravitational energy. The parameter f3 tells us whether clouds are 
gravitationally bound (/9 < 2), and whether in virial equilibrium (/9 ~ 1). The figure shows 
that the virialized clouds with Mc ^ 6 x 10^ M© have formed in model S1024. But in model 
R1024 the identified objects have Mc < 10^ M0, and none are gravitationally bound. We 
note that even in model RH512 with a half rotation speed, no gravitationally bound clouds 
were found. 

The bottom row of Figure 3 shows the specific angular momentum of identified clouds. 
At t = 4tcooi the specific angular momentum is different in the two models, although other 
properties are similar. As mentioned earlier, clumps formed in the early stage gain angular 
momentum through the Coriolis effect, so jc is higher in model R1024 than in model S1024. 
But during later stage angular momentum is acquired efficiently through torque and merging 
in model S1024, while it is not in model R1024. Hence for a given mass, jc becomes higher 
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in model S1024 than in model R1024. 



4. Summary 

We study the effects of rotation in protogalactic disk regions on the formation of 
structures via TGI in the protogalactic environment. A simplified setting was considered, 
where a gas of primordial composition evolves from initial density perturbations in a uni- 
formly rotating box. In Paper 1, we found that without rotation, virialized clouds of mass 
Mc ^ 6 X 10^ M0 can form as a result of TGI although they form with large angular momen- 
tum of spin parameter (Ag) ~ 0.3. In this paper we find that with rotation whose angular 
speed is comparable to that of the Galactic rotation at the Solar circle, Vto = 27 km s~^ kpc~ , 
the Coriolis force suppresses gravitational infall and merging and disperses the gas to fila- 
mentary and sheet-like structures. As a result, instead of massive virialized clouds formed 
in non-rotating models, clumps with M^ ^10^ M©, which are gravitationally unbound and 
often transient, are only found. 

We conclude that the rotation in protogalactic disk regions has destructive effects on 
the formation of clouds, rather than alleviate the angular momentum problem discussed in 
Paper 1. The results in this paper and those in Paper 1 combined imply that it would be 
difficult for globular clusters to have formed via TGI in protogalaxies, and even less likely 
in rotating disk regions. 

A few notes on our results: 1) Primordial gas is considered in this paper. But adding 
metal of order of 0.1 Zq, which would enhance thermal processes, does not change the results, 
because it is the gravitational processes that are responsible for the formation of massive 
clouds but suppressed by the Coriolis force. We confirmed it with another simulation with 
metalicity of 0.1 Zq, although we do not report the simulation in the paper. We point, 
however, that changing other parameters, for instance increasing the gas density by a factor 
of 10 in particular, would make differences. 2) The size of the simulation box we used, 
L = 10 kpc, along with periodic boundaries is too large to represent a protogalactic disk 
region. In order to check the effects of box size, we performed simulations in a smaller box 
of 2.5 kpc with 256'^ cells with and without rotation. Again although we do not report 
the simulation in details, the results confirmed that the formation of gravitationally bound 
objects is prohibited by the Coriolis force, regardless of the box size. 3) To be more realistic, 
differential rotation instead of uniform rotation needs to be investigated. 
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Table 1. Model Parameters for Simulations 



Model 


No. of Grid Zones 


^cnd (.^coolj 


fib 


S1024 
R1024 
RH512 


10243 
10243 
512^ 


16 
16 
20 


0.0 
O.Sfio 


^cool ~ 


= 2 X 10^ yrs 

27 km s"^ kpc"^. 
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Fig. 1. — Evolution of the density power spectrum in models S1024, R1024, and RH512. 
Circles represent the initial power spectrum at t = 0. Lines show the power spectrum at 
2 tcooi, 4 tcooi, 6 tcooi, • • •, 16 tcooi in models S1024 and R1024, and at 4 tcooi, 6 tcooi, • • •, 
20 tcooi in model RH512. 
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Fig. 2. — Isodensity surfaces inside a 2.5 kpc box of 256^ grid zones, which is (1/4)^ of the 
total simulation box, at 8 tcooi, 12 tcooi, and 16 tcooi in model S1024 (top panels), and in 
model R1024 (bottom panels). Green surfaces correspond to lOpo, yellow surfaces to lO^po; 
and red surfaces to lO^po- Here po is the mean initial density. The box is oriented in such 
a way that the x, y, and z— axes are from near to far, from bottom to top, and from left to 
right, respectively. Rotation is along the z— axes in model R1024. 
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Fig. 3. — Differential number of clouds dNc/d{\ogMc), energy ratio, 2{Et + Ek)/\Eg\, and 
specific angular momentum jc-, as a function of cloud mass M^., at four different times in 
models S1024 (blue color) and R1024 (red color). Here jc is in cgs units. 



